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Abstract: The OPERA experiment, exposed to the CERN to Gran Sasso beam, 
collected data from 2008 to 2012. Four oscillated Vr Charged Current interaction candidates 
have been detected in appearance mode, which are consistent with —)• Vr oscillations 
at the atmospheric Am^ within the “standard” three-neutrino framework. In this paper, 
the OPERA Vr appearance results are used to derive limits on the mixing parameters of a 
massive sterile neutrino. 
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1 Introduction 

The OPERA experiment [1] operated in the CERN Neutrinos to Gran Sasso (CNGS) 
beam produced at CERN and directed towards the Gran Sasso Underground Laboratory 
of INFN (LNGS), 730 km away, where the detector is located. The experiment is unique 
in its capability to observe appearance on an event-by-event basis. Nuclear emulsion 
films instrumenting the target allow the detection of the short-lived r lepton decay, and 
hence the identification of Ur Charged Current (CC) interactions. The standard three- 
neutrino oscillation framework predicts —)• oscillations with close-to-maximal mixing 
at the so-called atmospheric scale, Am |2 ~ 2.4 x 10“^ eV^ [2], i.e. in the oscillation 
parameters region discovered by detecting atmospheric neutrinos [3]. OPERA has observed 
four Vr CC interaction candidate events [4] [5] [6] [7], consistent with the expectation of the 
standard oscillation framework at this scale. This result represents the first direct evidence 
of Vr oscillation in appearance mode. 

In the present paper, limits are derived on the existence of a massive sterile neutrino. 
The excess of Ve (Pg) observed by the LSND [8] and MiniBooNE [9] collaborations and the 
so-called reactor [10] and Gallium [11] neutrino anomalies are also interpreted as due to the 
existence of a fourth sterile neutrino with mass at the eV scale. In relation to this issue, it 
is worth mentioning that the effective number of neutrino-like species decoupled from the 
primeval plasma measured by the Planck collaboration is 3.15 ± 0.23 at 95% Confidence 
Level (CL) [12]. 

Neutrino oscillations at large Am^ have been searched for by several short baseline 
experiments. The most stringent limits on —>■ Vr oscillations were set by the NOMAD 
[13] and CHORUS [14] experiments, with high sensitivity for Am^ values larger than 10 
eV2. 

In the following, a short description of the OPERA experimental setup and of the pro¬ 
cedure used to detect Ut interactions is given, the data analysis is described and exclusion 
regions in the parameter space are derived. 
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2 Detector, beam, and data sample 


In OPERA, CNGS neutrinos interacted in a massive target made of lead plates interspaced 
with nuclear emulsion films acting as high accuracy tracking devices [15]. This kind of 
detector is historically called Emulsion Cloud Chamber (ECC). The OPERA detector is 
made of two identical Super Modules, each consisting of a target section and a muon 
magnetic spectrometer. Each target section has a mass of about 625 tons and is made of 
~75000 units called bricks. A target brick consists of 56 lead plates Imm-thick interleaved 
with 57 nuclear emulsion films 300 /rm-thick for a mass of 8.3 kg. Its total thickness along 
the beam direction corresponds to about 10 radiation lengths. The bricks are assembled 
in vertical walls instrumented with scintillator strips {Target Tracker detectors, TT) to 
trigger the read-out and locate neutrino interactions within the target. 



Figure 1. Fluxes of the different CNGS beam neutrino components at LNGS. 


OPERA was exposed to the CNGS beam generated by protons from the SPS ac¬ 
celerator at CERN [16]. The contaminations of i^e and Pe CC interactions at LNGS, 
relative to CC interactions, are 2.1%, 0.9% and less than 0.1%, respectively. The con¬ 
tamination of prompt i^r is negligible. The average energy is 17 GeV. The energy spectra 
of the four beam components at the detector site are shown in figure 1 [17]. 

The data taking, based on a minimum bias interaction trigger from the TT scintillators, 
started in 2008 and ended in December 2012. OPERA collected data corresponding to 
17.97 X 10^® protons on target (pot) with 19505 recorded events. The data sample used in 
this analysis is defined following the selection criteria described in [5] and corresponds to 
about 75% of the total statistics. 
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3 Search for Vr interactions 

Bricks selected as candidates to contain CNGS neutrino interactions are analysed following 
the procedure described in detail in [5]. Here we just recall the main steps of the analysis. 

The brick where a neutrino interaction occurred is predicted by the electronic detectors 
and extracted from the target by an automatic brick manipulator system. Two extra low 
background emulsion films {Changeable Sheets, CS) [18] located downstream of the brick 
act as an interface between the brick and the electronic detectors. If the measurement of 
the CS yields tracks related to the neutrino interaction, the emulsion films of the brick are 
developed. Their analysis provides the three dimensional reconstruction of the neutrino 
interaction and of the possible decay vertices of short-lived particles [19] with micrometric 
accuracy. 

This procedure has led to the detection of four Uj- CC interaction candidates. The total 
expected background in the analysed sample amounts to 0.23 ± 0.05 events. The absence 
of a —>• Vt oscillation signal, i.e. the hypothesis of the four events being background, is 
excluded with a significance of 4.2 a [7]. 


4 Sterile neutrino search via —)• oscillations 

In [7] the detection of four CC events is compared to the expectation for —)• 

oscillations in the atmospheric sector, computed within a simplified two-flavour scheme 
assuming full mixing and jAm 32 l = 2.32 x 10“^ eV^ [20]. The expected number of events 
is 2.30 ± 0.46 (2.21 ± 0.44) assuming normal (inverted) hierarchy of neutrino masses; the 
number is obtained by rescaling the value given in [7] for jAm| 2 l ~ ~ [Am^j = 

2.43 X 10“^ eV^ (jAm^j = 2.38 x 10“^ eV^) [2], where Am^ is defined as _ 

By including the background, 2.53 ± 0.46 (2.44 ± 0.44) events are expected in total. The 
error, which is dominated by the uncertainty on the r detection efficiency and on the 
Vr interaction cross section, also takes into account the experimental precision on the 
atmospheric oscillation parameters. The observation of four events is compatible with 
these expectations. Despite the limited statistics, an excess or a deficit of Ur interactions 
due to — 7 - Vt oscillations induced by the mixing with a sterile neutrino can be evaluated. 

In presence of a fourth sterile neutrino with mass m 4 , the oscillation probability is a 
function of the 4x4 mixing matrix U and of the three squared mass differences. Defining 
C = 2\U^'i\\Ur‘i\, ^ij = 1-27 Anifj LjE {i,j = 1,2,3,4), = Arg{Uf, 3 U*^U*^Ur 4 ) and 

sin20^T- = 2\Uf^4\\Uri\, the —>■ Ur oscillation probability P(E) can be parametrised as: 


P{E) = sin^ A 31 -|- siT?29^r sin^ A 41 

- C sin 20^,- cos(f>nr sin 2 A 3 i sin 2 A 4 i 

— C sin 26HT sin (1)^^ sin^ A 31 sin 2 A 41 
-|- 2 C sin 29HT cos sin^ A 31 sin^ A 41 
-|- C sin 29 sin sin 2 A 31 sin^ A 41 


(4.1) 
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where and /^m\i are expressed in eV^, L in km and E in GeV. Given the long baseline 

and the average CNGS neutrino energy, P(E) is independent of Am^i, since A 21 ~ 4x 10“^. 
The terms proportional to are CP-violating, while those proportional to sin 2 A 3 i 

are sensitive to the mass hierarchy of the three standard neutrinos, normal (Aml^ > 0 ) or 
inverted (Am|^ < 0). Matter effects have been checked to be negligible for Am|^ > 1 eV^. 

Observed neutrino oscillation anomalies [ 21 ], if interpreted in terms of one additional 
sterile neutrino, suggest |Am|;^| values at the eV^ scale (the so-called 3-1-1 model). In the 
following, unless stated otherwise, the analysis will be restricted only to positive Aml^ 
values, since negative values are disfavoured by results on the sum of neutrino masses 
from cosmological surveys [12]. For Am^^^ > 1 eV^, at the concerned domain of L/E and 
taking into account the finite energy resolution, sin 2 A 41 and sin^ A 41 average to 0 and 
respectively. The oscillation probability P{E) can thus be approximated to [21]: 


P{E) = sin^Aai -h ^ sin^ 20^^ 

-|- C sm.26^T- cos (j)^r sin^ A 31 (4.2) 

-I- ^ C sin 20^^ sin 4)nr sin2A3i. 

In order to obtain an upper limit on sin^ 26 ^r at high values of Am^^, the likelihood 
is defined as L{(j)nT, sin^ 20^,-, G^)= [Eln\, where n = 4 is the number of Ur candidate 
events and p is the expected number of events, fj, = rib + A f 4>{P)P{P)^{P)^[P) dE. 
rib = 0.23 is the expected number of background events [7], 4){E) is the flux shown in 
figure 1, P{E) is the oscillation probability given by equations (4.1) or (4.2), cr{E) is the 
Ur GG interaction cross section, €{E) is the r detection efficiency and A is a normalisation 
factor proportional to the fraction of the analysed sample and to the target mass. 

The analysis presented here is based on the asymptotic distribution of the log 

likelihood ratio test statistics: q = —2 ln(T((^^T-, sin^ 20^T-)/Lt)), where Lg = ri^ln\ 
and L{(j)nr, siv? 29 ^t) is the profile likelihood obtained by maximising L{(pnT ,sin^ 20 ^,-, G^) 
over By definition, ranges between 0 and 1, but for any pair of values of sin^ 26 
and (piiT, it is limited by the unitarity of the mixing matrix; the likelihood is maximised 
accordingly. The value of jAm| 4 j has been fixed to 2.43 x 10“^ eV^ for the normal hierarchy 
and to 2.38 x 10“^ eV^ for the inverted hierarchy of the three standard neutrinos [2]. 

In figure 2(a) the 90% CL exclusion limits are presented for both normal and inverted 
mass hierarchies in the parameter space of (jr^T vs sin^ 29 The edge of the excluded region 
ranges from 0.088 to 0.136 in sin^ 20^,- for both mass hierarchies of the three standard 
neutrinos. For any fixed value of Q is distributed according to a statistics with one 
degree of freedom. Profiling the likelihood also over (pitr, as shown in figure 2(b), an upper 
limit of 0.116 is obtained at 90% CL on sin^ 20^,-, almost independently of the hierarchy of 
the three standard neutrino masses. A negligible difference arises from the different j 
value used in the analysis. The sin^ 29 upper limit is affected by a 20% systematic error 
from the uncertainties on the r detection efficiency and Ur interaction cross section. 
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Figure 2. (a) 90% CL exclusion limits in the vs sin^ parameter space for normal (NH, 

dashed red) and inverted (IH, solid blue) hierarchies assuming > 1 eV^. Bands are drawn 

to indicate the exclnded regions, (b) Log likelihood ratio as a function of sin^ 20^,- for = 0 
(dashed line) and for the profile likelihood (continuous line). 



Figure 3. 90% CL exclusion limits (blue line) in the |C/i- 4 p vs plane assuming Am\^ > 

1 eV^. The unitarity bound (grey line) is also shown. Bands are drawn to indicate the excluded 
regions. 


Given the definition of sin^ in terms of C /^4 and Ur a, it is possible to translate the 
upper limit on sin^ 20 ^,- into an exclusion curve in the vs |C/t- 4 P plane, as shown in 

figure 3 together with the unitarity bound {\U^a\^ + |C/r 4 p < !)• 

To extend the search for a possible fourth sterile neutrino down to small values, 
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the likelihood has been computed using the GLoBES software [23], in order to take into 
account also matter effects and the non-zero value of Amli- The likelihood has been 
profiled also on the value. More details on the analysis are available in [22], In 

figure 4 the 90% CL exclusion plot is reported in the Aml^ vs sin^ 20^,- parameter space. 
The most stringent limits of direct searches for —)• Ur oscillations at short-baselines 
obtained by the NOMAD [13] and CHORUS [14] experiments are also shown. Our analysis 
stretches the limits on Am|^ down to 10 “^ eV^, extending the values explored with the r 
appearance searches by about two orders of magnitude at large mixing, for sin^ 26 ^ 0.5. 
For maximal mixing, the 90% CL excluded region extends down to Amh = 7.4 (5.2) 
xl 0 “^ eV^ for normal (inverted) hierarchy of the three standard neutrino masses, with a 
10 % systematic error deriving from the uncertainties on the r detection efficiency and Ur 
interaction cross section. 
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Figure 4. OPERA 90% CL exclusion limits in the vs sin^ 20^^ parameter space for the 

normal (NH, dashed red) and inverted (IH, solid blue) hierarchy of the three standard neutrino 
masses. The exclusion plots by NOMAD [13] and CHORUS [14] are also shown. Bands are drawn 
to indicate the excluded regions. 

A narrow region is excluded at 90% CL at Amh ~ 10“^ eV^ for the normal hierarchy 
of the three standard neutrinos. It arises from a suppression of the —)• Vr oscillation 
probability due to the presence of the sterile neutrino. Instead, the oscillation probability 
is enhanced at full mixing and high Aml^ values. For a number of r neutrino candidates 
equal to the expectation in the three neutrino framework, the excluded region at Am|g ~ 
10 “^ eV^ would disappear. 

The analysis was performed assuming Am|^ > 0. Since present limits on the sum of 
neutrino masses from cosmological surveys do not exclude small negative values for Am\i, 
the analysis was repeated following this assumption. The exclusion plots obtained in this 
way are similar to those of figure 4, but with hierarchies exchanged. It is worth underlining 
that the results obtained in the 3-1-1 model, shown in figures 2 and 3, are independent of 
the sign of Amh, as is the probability in equation 4.2. 

Assuming CP conservation, that implies sin = 0, and neglecting terms in sin^ A 31 
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in equation (4.1), of the order of 10“^ at CNGS energies Ecngs, the oscillation probability 
at high values of approximates to that of a two-flavour model parametrised in terms 

of two effective mixing parameters, 0e// and 

lim P{E\CP conservation) ~ sin^ ' sm^(1.27 Avniff L/E) 

E^Ecngs 

In this framework, with 4 observed events and 2.53 ± 0.46 events expected from the normal 
hierarchy of standard oscillations, including background, the upper limit on the number of 
additional i^r events, evaluated in the Feldman-Cousins approach [24], is 6.4 at 90% CL. 
The upper limit on sin^ “^^eff is 0.069, to be compared with 0.076 for sin^ 20 ^t- at (jy^r = 0 
(see hgure 2(b)). 

5 Conclusions 

The OPERA experiment was designed to observe —>• Vr oscillations through I'r ap¬ 
pearance at a baseline of 730 km in the CNGS beam. Exploiting its unique capability to 
identify r neutrino interactions, OPERA has observed four I'r CC candidate interactions, 
consistent with the expected number of oscillation events in the standard three-neutrino 
framework, 2.30 ± 0.46 (2.21 ± 0.44), for the normal (inverted) mass hierarchy and 0.23 
± 0.05 background events. 

In this paper we present limits on the existence of a sterile neutrino in the 3-1-1 neutrino 
model. At high values of Am|^, the measured 90% CL upper limit on the mixing term 
sin^ 29 = 4|t/^4p|t/T-4p is 0.116, independently of the mass hierarchy of the three standard 
neutrinos. The OPERA experiment extends the exclusion limits on Am\i in the —)> Vr 
appearance channel down to values of 10“^ eV^ at large mixing for sin^ 20^,- > 0.5. 
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